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ABSTRACT 

Aims. We study the influence of large-scale photospheric motions on the destabilization of an eruptive filament, observed on October 
6, 7, and 8, 2004, as part of an international observing campaign (JOP 178). 

Methods. Large-scale horizontal flows were invetigated from a series of MDI full-disc Dopplergrams and magnetograms. From the 
Dopplergrams, we tracked supergranular flow patterns using the local correlation tracking (LCT) technique. We used both LCT and 
manual tracking of isolated magnetic elements to obtain horizontal velocities from magnetograms. 

Results. We find that the measured flow fields obtained by the different methods are well-correlated on large scales. The topology of 
the flow field changed significantly during the filament eruptive phase, suggesting a possible coupling between the surface flow field 
and the coronal magnetic field. We measured an increase in the shear below the point where the eruption starts and a decrease in shear 
after the eruption. We find a pattern in the large-scale horizontal flows at the solar surface that interact with differential rotation. 
Conclusions. We conclude that there is probably a link between changes in surface flow and the disappearance of the eruptive 
filament. 
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1. Introduction 

Dynamic processes on the Sun are linked to the evolution of the 
magnetic field as it is influenced by the different layers from 
the convection zone to the solar atmosphere. In the photosphere, 
magnetic fields are subject to diffusion due to supergranular 
flows and to the large-scale motions of differential rotation and 
meridional circulation. The action of these surface motions on 
magnetic fields plays an important role in the formation of large- 
scale filaments (Mackay and Gaizauskas 2003). In particular, the 
magnetic fields that are transported across the solar surface can 
be sheared by dynamic surface motions, which in turn result in 
shearing of the coronal field. This corresponds to the formation 
of coronal flux ropes in models, which can be compared with 
Ha filament observations (Mackay and van Ballegooijen 2006b). 
Many theoretical models try to reproduce the basic structure and 
the stability of filaments by taking surface motions into account, 
as quoted above. These models predict that magnetic flux ropes 
involved in solar filament formation may be stable for many 
days and then suddenly become unstable, resulting in filament 
eruption. Observations show that twisting motions are a very 
common characteristic of eruptive prominences (see for exam- 
ple Patsourakos and Vial 2002). However, it is still unknown 
whether the magnetic flux ropes emerge already twisted or if it is 
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only the photospheric motion that drives the twisting of the fila- 
ment magnetic field. The destabilization of the filament can also 
be linked to oscillations (Pouget 2006). Therefore, it remains un- 
certain as to the mechanisms that drive filament disappearance. 
Destabilization can come from the interior of the structure or by 
means of an outside flare. 

In a previous paper (Rondi et al. 2007, hence forth Paper I), 
local horizontal photospheric flows were measured at high spa- 
tial resolution (0.5") in the vicinity of and beneath a filament be- 
fore and during the filament's eruptive phases (the international 
JOP178 campaign). It was shown that the disappearance of the 
filament originates in a filament gap. Both parasitic and normal 
magnetic polarities were continuously swept into the gap by the 
diverging supergranular flow. We also observed the interaction 
of opposite polarities in the same region, which could be a can- 
didate for initiating the destabilization of the filament by causing 
a reorganization of the magnetic field. 

In this paper we investigate the large-scale photospheric 
flows at moderate spatial resolution (2") beneath and in the 
vicinity of the same eruptive filament. The observation and 
coalignment between data from various instruments are ex- 
plained in Section 2. In Section 3, we describe the different meth- 
ods of determining the flow field on the Sun surface. The large- 
scale flows associated with the filament are shown Section 4. 
The properties of these flows before and after the filament erup- 
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Fig. 1. MDI longitudinal magnetic field data on October 7, 2004, 
at the beginning of the sequence 9:40 UT and at its end 22.40 UT. 
In order to see the magnetic field evolution, we overplot a con- 
tour of filament observed at 13:30 UT in both figures. 



tion are described in Section 5. In Section 6, we investigate the 
topology of horizontal flows in the filament area over the 3 days 
around the filament eruption. A discussion of the results and gen- 
eral conclusions can be found in Section 7. 



2. Observations 

During three consecutive days of the JOP 178 campaign, Oct 6, 
7, and 8, 2004 (http://gaia.bagn.obs-mip.fr/jopl78/index.html), 
we observed the evolution of a filament that was close to the cen- 
tral meridian. We also observed the photospheric flows directly 
below the filament and in its immediate vicinity. The filament ex- 
tends from -5° to -30° in latitude. A filament eruption was ob- 
served on October 7, 2004, at 16:30 UT at multiple wavelengths 
from ground and space instruments. The eruption produced a 
coronal mass ejection (CME) at approximately 19:00 UT that 
was observed with LASCO-2/SOHO and two ribbon flares ob- 
served with the SOHO/EIT. MDI/SOHO longitudinal magnetic 
field and Doppler velocity were recorded with a cadence of one 
minute during the 3 days (see Table 1). The Air Force O-SPAN 
telescope located at the National Solar Observatory/Sacramento 
Peak provided a full-disc Ha image every minute. The pixel 
sizes were 1.96" for MDI magnetograms and Dopplergrams and 
1.077" for O-SPAN Ha images. Table 1 summarises the charac- 
teristics of all the observations of JOP 178 used in our analysis. 

We coaligned all the data obtained by the different in- 
struments (see Paper I for a complete description of the co- 
alignment procedure). Our primary goal was to derive the hor- 
izontal flow field below and around the filament. Co-alignment 
between SOHO/MDI magnetograms and O-SPAN data was ac- 
complished by adjusting the chromospheric network visible in 
Ha (O-SPAN) and the amplitude of longitudinal MDI magne- 
tograms to an accuracy of one pixel (1.96"). 

The general magnetic context before and after the filament 
eruption is shown Fig.Q] 



3. Determination of the photospheric flows 

Horizontal flows on the solar surface may be measured through 
the proper motion of the plasma or by the effects of the plasma 
on the magnetic structures. 



3. 1 . Dopplergram processing 

In order to map the horizontal component of the large-scale 
photospheric plasma velocity fields, we applied local corre- 
lation tracking (LCT; November 1986) to a set of full-disc 
Dopplergrams obtained by the MDI instrument onboard SoHO. 
The aim of this method is to track the proper motion of su- 
pergranules that are clearly detectable on Dopplergrams every- 
where except for the disk centre. The Dopplergrams were pro- 
cessed following the procedure described in Svanda et al. (2006) 
with slight modifications for our data. Hereafter we refer to this 
method as LCT-Doppler. Initially we suppress the solar p-mode 
oscillations. Using a weighted temporal average (see Hathaway 
1988) described by the formula: 
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where Af is the time interval between a given frame and the cen- 
tral one (in minutes), b — 16 min and a — 8 min. The normal- 
ized version of this function has been applied to the data. This 
filter reduces the amplitudes of the solar oscillations in the 2- 
4 mHz frequency band by a factor of more than five hundred. 
The oscillations in each frame were reduced using a window 
of 31 successive frames. The different time series were tracked 
using the Carrington rotation rate (with an angular velocity of 
13.2 degrees per day), so that all the frames have the same he- 
liographic longitude of the central meridian (Iq = 62.24 °). The 
tracked data were remapped into a sinusoidal pseudocylindrical 
coordinate system (also know as Sanson-Flamsteed grid) to re- 
duce the distortion of structures in the Dopplergrams caused by 
the geometrical projection to the disc. The sinusoidal projection 
is suitable to describe the behaviour on the large scales. Tracked 
and remapped time series then undergo a k-a> filtering with cut- 
off velocity of 1500 m s _1 to suppress the noise coming from the 
groups of granules and the change of contrast of supergranular 
structures due to the solar rotation. The individual frames were 
apodizedby 10% using a smooth function, the same apodization 
took place in the temporal domain. The resulting data series of 
tracked, remapped and filtered frames were then ready for track- 
ing. 

The LCT method applied to full-disc Dopplergrams is char- 
acterised by a Gaussian correlation window (FWHM = 60") 
and a time lag between correlated frames of 1 hour (basically 60 
frames). In all cases, one half of the intervals were before the 
eruption and the second half after the eruption. All the pairs of 
correlated frames in the studied intervals were averaged to in- 
crease the signal to numerical noise ratio. 

3.2. Magnetogram processing 

The second method by which we determine motions on the 
solar surface used the full-disc magnetograms obtained by 
MDI/SoHO. To reduce the distortion of structures seen in the 
magnetograms caused by the geometrical projection to the disc, 
we applied Sanson-Flamsteed grid projection. To measure the 
differential motions of features on the solar surface, the data 
were aligned on a band along the equator. Due to the numer- 
ous magnetic structures to be tracked manually, the field of view 
was limited to -35° to 20° in longitude and -3° to -30° in lat- 
itude. The displacement of the longitudinal magnetic structures 
visible on the magnetograms (both positive and negative) was 
determined with two different methods: 
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Telescope 



Datatype 



Field of view Pixel size Cadence 



Time U.T 



1SOON 



Ha 



full-disc 



1.077" 



1 min 



14:05- 
13:30- 
13:37- 



22:35 October 6 
20:31 October 7 
22:35 October 8 



MD1/SOHO magnetogram full-disc 1.96" 1 min 20:49-23:49 October 6 

Doppler velocity 9:44-22:50 October 7 

1 min 6:56-12:52 October 8 

1 min 16:09-19:33 October 8 



1 . The first approach, hereafter named manu-B, was to manu- 
ally locate each magnetic structure in each magnetogram to 
determine its trajectory and its horizontal velocity. Once the 
velocities in the field of view were determined, as the mag- 
netic structures do not sample the field of view unifomly, 
we applied a reconstruction of the velocity field based on 
multi -resolution analysis described by Rieutord et al. (2007), 
which allows us to limit the effects of the noise and error 
propagation. 

2. The second approach, named LCT-B, was to apply the LCT 
on magnetic structures with absolute values greater than 
25 Gauss to reduce the noise. The LCT method applied to 
these magnetograms used a Gaussian correlation window 
(FWHM = 60") and a time lag between correlated frames 
of 1 hour. 

The horizontal flow field measured using the plasma 
(Dopplergrams) differs from that measured using magnetic fea- 
tures. This results because the magnetic field structures are not 
actually passive scalars; they can interact with the plasma and are 
also constrained by their interactions in the upper atmosphere. It 
is therefore not surprising to observe small differences in the ve- 
locity fields derived from the various tracers (Dopplergram or 
magnetic structures). 

4. Photospheric flow pattern below and around the 
filament 

In this Section we describe the flows associated with the filament 
eruption, with particular emphasis on the filament evolution and 
properties of the the mean East- West (zonal) velocities. 

4.1. Flow fields 

The flow fields in the vicinity of the filament obtained from 
13-hour averages of velocities found using the three different 
methods described above are shown in Fig. [2] Inspection of 
these maps shows that all three methods provide similar large- 
scale velocity patterns. The amplitudes of the flows are given in 
Table 2. The LCT method applied to magnetic structures with 
amplitude greater than 25 Gauss appears to underestimate the 
amplitude of the flows obtained from the two other methods. 
This is partially because of to the large window (FWHM = 60") 
used in the LCT method, combined with the uneven distribution 
of the magnetic structures in the field of view. The correlation 
between zonal component, v x , between the three methods is be- 
tween 0.48 and 0.40, while the correlation for meridian com- 
ponent, v y , is between 0.21 and 0.11. The low correlation co- 
efficients for the meridian component probably occurs because 
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Fig. 2. The horizontal photospheric flow fields that were derived 
using the various methods described in the text. Upper panel - 
manu-B method; middle panel - LCT-Doppler method; and bot- 
tom panel - LCT-B method. The filament as observed by ISOON 
in H-halpha on October 7, 2004, is superimposed. The hexagon 
indicates the location where the filament eruption started. 
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the main direction of the flows in our field of view is zonal and 
theamplitude of the meridian component is small. We estimate 
that the directional error in our measurements of the velocities 
is about ±5°. A small error in determining the direction of the 
flow, for example 10°, can affect the amplitude of the meridian 
component by a factor two. This error would greatly decrease the 
correlation in the meridian component of the flow. However, as 
seen in Fig. [2] the general trend is similar in the velocity fields re- 
sulting from the three different methods. In particular, the north- 
south stream that disturbs the differential rotation around -25° is 
easily visible, and most of the large-scale features of the vector 
orientations can be identified. We note that the lowest correlation 
coefficient is found between LCT-B and the other methods manu- 
B, LCT-Doppler. This results when the LCT-B method clearly 
underestimates the amplitude flows. In agreement with the pre- 
vious results of Schuck (2006), we conclude that this method is 
probably not suitable for accurately estimating horizontal veloc- 
ities of magnetic footpoints. 

The correlation is higher in the longitudinal region be- 
tween heliographic coordinates of -5° and 20° where it is 
between 0.58 and 0.41 for v x and between 0.34 and 0.17 for 
v y . In this region, the large-scale flows are well-structured and 
show both converging and diverging velocity patterns. We 
observe in particular a large-scale stream in the north-south 
direction parallel to the filament located about 10° to the 
east and between -20° and -30° in latitude and 58°and 47° 
in longitude. This flow stream is clearly visible in Fig. [2] 
(upper and middle panels), and its dynamics can be seen 
at http://gaia.bagn.obs-mip.fr/jopl78/oct7/mdi/7oct-mdi.htm 
Around -20° in latitude, the velocities of the differential rotation 
amplitude appears to dominate. However, in both measurements 
we observe that the north-south large-scale stream on the 
eastern edge of the filament disturbs the regular differential 
rotation. The north-south stream is located close to where the 
filament eruption begins (longitude I = 56°, latitude b = -26° 
in Carrington coordinates). In the manu-B and LCT-B methods, 
the stream appears closer to the filament. The amplitudes of 
the southward motions are manu-B 31.2 ms 1 , LCT-Doppler 
40 ms~', and LCT-B 13 ms~' which are close to the mean 
observed flows (see Table 2). Below the latitude of -20°, 
the combination of differential rotation and the north-south 
stream cause opposite polarities to move closer together, which 
strongly increases the tension in the magnetic field very close 
to the starting point of the filament eruption. Our measurements 
show a good agreement between the manu-B and LCT-Doppler 
methods. 

4.2. Filament evolution 

The filament's evolution can be seen in Fig. [3] The north-south 
stream flow visible in Fig. [6] (left) crosses over the part of the 
filament labeled A. The arrow on Fig.[3]indicates the same fixed 
point (325", 167") in all of the subframes. We observe a gen- 
eral southward motion of both the A and B segments of the fil- 
ament. More precisely, we measure a tilt of these two filament 
segments at the point of their separation. Between 16:07 UT and 
16:58 UT, the long axis of segment A rotates by an angle of 
12° (clockwise) relative to its western end, and the long axis of 
segment B of the filament rotates by an angle of 5.5° (clock- 
wise) relative to its western end. These rotations are compatible 
with the surface flow shown in Fig. [6] (left) and in particular the 
north-south stream flow. We did not find a singular pivot point 
where differential rotation did not displace the filament with re- 
spect to the flare location (Mouradian 2007) for the present fila- 
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Fig. 3. Evolution of the filament during its eruption around 
16:30 UT on October 7, 2004. The arrow indicates a fixed po- 
sition for all the subframes. A and B denote two parts of the 
filament. 



Table 2. Velocity amplitude in m s 





manu-B 


LCT-Doppler 


LCT-B 


mean velocity 


43.6 


28.7 


10.6 


maximum velocity 


111.1 


87.0 


32.5 


minimum velocity 


0.9 


0.26 


0.04 



ment. The southern segment of the filament does not reform after 
the eruption. This sudden disappearance shows there was an im- 
portant change in the sun and not only in the solar atmosphere 
(Mouradian 2007). 

4.3. Mean zonal velocities 

Figure |4] show plots of the mean zonal velocities resulting from 
the three different velocity determination methods as functions 
of latitude. The mean zonal velocities from the manu-B method 
clearly show the differential rotation profile with a plateau 
around -25°, which corresponds to the effects of the north-south 
stream discussed above. A similar profile, but with a lower am- 
plitude, is visible in the mean zonal velocities obtained from the 
LCT-B method. These two methods measure the displacement of 
the magnetic structures on the Sun's surface. The LCT-Doppler 
method, which measures the motion of the photospheric plasma, 
shows a mean zonal velocity profile in which differential rota- 
tion is clearly visible, along with a strong secondary maximum 
visible at -23° of latitude. The secondary maximum indicates 
a decrease in the amplitude of the v x component because the 
flows in that region are oriented more in the north-south direc- 
tion. That is partly due to the presence of the north-south stream 
described above and to the local organisation of the flow. In par- 
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Fig. 4. Mean zonal velocity is computed between -35° and +20° 
in longitude and -30° to 0° in latitude in the Carrington system. 
Dash-dot line and diamonds are zonal velocity for the motions 
measured by hand on B long maps (manu-B method). Solid line 
and squares represents are zonal velocity for the LCT applied on 
Doppler images {LCT-Doppler method). The dashed line and 
triangles represents mean zonal velocity for the LCT method 
applied on B long maps (LCT-B method). 



u 

C: 
> 
O 

c 
o 



cz 

D 
Cl) 




-30 -20 -10 

Heliographic latitude [deg] 

Fig. 5. Profiles for 7 Oct in different longitudinal belts solid line 
and squares -25° to -17° degrees with respect to the centre of 
the disc (not Carrington coordinates), dashed line and triangles 
for 0° to 20°. 



ticular, converging and diverging flows in this region seem to 
have a North-South orientation. 

To distinguish the effects of the north-south stream from dif- 
ferential rotation, we computed the mean zonal velocities from 
-25° to -17° in longitude centred on the region where the north- 
south stream is present and from 0° to 20° in longitude where, 
differential rotation dominates, with respect to the disc centre ve- 
locity. The mean zonal velocities in the longitudinal belt, where 
the north-south stream is visible, clearly exhibits a secondary 
maximum (Fig. [5]l indicating that the solar rotation rate at this 
location is closer to that of the equator. The mean zonal ve- 
locities computed in the standard differential rotation belt show 
the classic latitude profile with a constant decrease down to the 
low latitudes. As a consequence, the plasma in the north-south 



stream, which transports magnetic structures, rotates faster at 
about -23°latitude than do the magnetic structure located in the 
belt of longitude between 0° to 20°. The combination of these 
different surface motions (stream and differential rotation) tends 
to bring together fields with opposite polarities, and this in turn 
constrains the magnetic field lines. 

We noted that the location of the starting point of the filament 
eruption is around -26° in latitude, which is very close to the 
secondary maximum in the mean zonal velocity. Thus surface 
motions that bring opposite polarities together may play a role 
in triggering the filament eruption. 



5. Flow fields before and after the eruption 

In this section we discuss the properties of the flow field just 
before and after the filament eruption, which occurred at about 
16:30 UT. Due to the length of the sequence and because it is 
easier to accurately track the flow using Doppler images than 
it is to track the small number of magnetic features above 25 
Gauss, only the measurements obtained with the LCT-Doppler 
are used in this section. 

At the point where the filament eruption begins (/ = 56°, 
b = -26° in Carrington coordinates), we detected a steepen- 
ing of the gradient in the differential rotation curve. During the 
eruption, the gradient flattens out and a dip forms. While differ- 
ential rotation curves describe mean zonal velocities over most 
of the disc, their change in gradient signifies a change in the 
stretching forces influencing the magnetic field loops over the 
area under study. We can express the surface rotation as an 
even power of sin <p : R - A + B sin 2 <p + C sin 4 <f>, where A 
is the angular velocity rate of the equatorial rotation and (f> the 
heliographic latitude. From the data we find that the constant 
values (with their errors in parentheses) are before the erup- 
tion A = 13.375(0.010), B = -1.46(0.10), C = -1.42(0.20) 
and after the eruption A = 13.404(0.010), B = -1.78(0.10), 
C = -1.24(0.20). All of the rates are synodical in deg day" 1 . 
The full-disc profiles did not change significantly from before to 
after the eruption. For example, for a latitude of -30° the zonal 
velocity has values of 12.92 (resp. 12.88) deg day" 1 (-34 ms" 1 , 
resp. -39 ms" 1 in the Carrington coordinate system), for a 
latitude of -20° the values are 13.18 (-2 ms" 1 ), resp. 13.18 
(-3 ms" 1 ), deg day" 1 . Although the parameters of the smooth 
fitted curve did not change too much, the local residual with 
respect to the smooth curve changed at the latitude where the 
filament eruption starts. 

Figure [6] displays the horizontal flows over a wide field of 
view measured using the LCT-Doppler method and then averag- 
ing the resulting velocities over 3 hours, before and after the fila- 
ment eruption. Before the eruption we can clearly see the north- 
south stream parallel to and about 10° east of the filament. This 
stream disturbs differential rotation and brings plasma and mag- 
netic structures to the south. Although differential rotation tends 
to spread the magnetic lines to the east, the observed north-south 
stream tends to shear the magnetic lines. After the eruption, only 
a northern segment of the filament is visible and the north-south 
stream has disappeared. 

To quantify the evolution of the horizontal flow before and 
after the filament eruption, we computed the change in the di- 
rection of the velocity vectors. The noise discussed in sect. 4.1, 
which can affect the direction of small magnitude vectors, tends 
to reduce the correlation between flows. In order to mitigate this 
error, we computed the magnitude weighted cosine of the direc- 
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Fig. 6. Horizontal motions measured by LCT-Doppler before (left) and after (right) the eruption over a wide field of view. The 
filament observed by ISOON on 2004 October 7 at 13:30 UT is superimposed 
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Fig. 7. Magnitude-weighted cosine map computed between flows before and after of the eruption of the filament. White areas 
indicate a good agreement between the fields before and after the filament eruption. Black areas indicate the location where a lot of 
changes between horizontal velocities take place. The filament observed by ISOON on 2004 October 7 at 1 3 : 30 UT is superimposed. 
The hexagon indicates the location where the filament eruption started. 



tion difference (as in Svanda et al., 2007), which is robust to the 
presence of noise. This quantity is given by the formula: 



where a and b are vector fields, a-b is a scalar multiplication and 
|a| is the magnitude. The closer this quantity is to 1, the better the 
alignment between two vector fields. 

Figure|7]displays the magnitude-weighted cosine map com- 
puted between flows before and after the eruption of the filament 
shown in Fig. [6] This map was computed by using a sliding win- 
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Fig. 9. The evolution of the velocity shear in zonal components 
in time. The eruption of filament took place at 16:30 UT 

dow with a size of 8.8"on a side (41.2"being the side of the data 
plot). The magnitude weighted cosine map (Fig. [7} reveals that 
changes in the vicinity of the north-south velocity stream are 
significant, while the horizontal flow field in the remainder of the 
field of view is more stable. Although variations in the flow field 
were expected to be more or less random over the field of view, 
we observe in that particular case that most of the variations be- 
tween before and after of the filament eruption are located in 
the north-south stream. The disappearance of the north-south 
stream after the eruption could be linked to the eruption or to a 
natural evolution of the photospheric flows. 

Figure [8] shows the flow field in more detail at the site where 
the eruption starts. The shear in the zonal component at the point 
where the eruption starts (/ = 56°, b - -26° in Carrington co- 
ordinates) is clearly visible and exits before and after the erup- 
tion, although the shape of the apparent vorticity has changed. 
This location corresponds to the area of upflow observed in the 
Meudon Ha Dopplergram (Fig. 8 in Paper I). 

We defined the shear as a difference between the mean zonal 
component v x in the area just North and just South of the point 
at which the eruption appeared to start. We obtained the mean 
zonal flow by averaging over boxes 2.3° on a side located 2.9° 
North and South of this point. The evolution in the shear veloc- 
ity computed as the difference between the mean flow in the two 
boxes as a function of time is shown in Fig. [9] Six 2-hour aver- 
ages of the flow fields were used to create this figure. The error 
bars come from estimating the noise measured on the basis of 
synthetic data from Svanda et al. (2006). One can see that the 
shear velocity increased before the eruption and decreased after 
the eruption. One hour before the eruption, the shear reached the 
value of (120+15) ms _1 over a distance of 5.2° (62000 km in 
the photosphere). After the filament eruption, we observed the 
restoration of ordinary differential rotation below 30° south. 

6. Evolution over 6, 7, and 8 October, 2004 

In Fig. [10] we compare the flow field in the filament region on 
the day of the eruption (October 7, 2004), with the flow fields 
on the preceding and following days. We see that the topology 
of the flows in this region changed over the three days: the daily 
evolution of the mean zonal profiles is shown in Fig. [TOl (bottom 
right). The dashed line with triangles is the mean zonal velocity 
for 6 October, 2004. This profile is relatively flat probably due to 
the short time sequence as only 3 hours of data were available. 
The differential rotation profile for October 7 shows a secondary 
maximum around -23° in latitude, as discussed in Sect. 4. The 
October 8 profile exhibits similar trend, but with a smaller am- 



plitude and an eastward velocity for latitudes greater than 10°. 
The secondary maximum appears strongly reduced, indicating 
restoration to a more regular differential rotation pattern in that 
zone. 

One day before the eruption, shear began to form at the site 
where the filament eruption is triggered (I = 56°, b = -26° in 
Carrington coordinates). The north-south stream is also visible. 
Both phenomena may store free energy in the coronal magnetic 
field configuration. The topology of the flow and the stream are 
different a day after the filament eruption suggesting that, after 
the disappearance of the southern part of the filament the con- 
ditions in the photosphere below, the filament became more re- 
laxed. This may suggest the mutual coupling of the photospheric 
flow and the configuration of the coronal magnetic field. To con- 
firm this idea, high-cadence high-resolution images and magne- 
tograms covering the eruption time would be needed. 

7. Discussion and conclusion 

Filaments or prominences are important complex structures of 
the solar atmosphere because they are linked to CMEs, which 
can influence the Earth's atmosphere and near-space environ- 
ment. Surface motions acting on pre-existing coronal fields play 
a critical role in the formation of filaments. They appear to 
reconfigure existing coronal fields, by twisting and stretching 
them, thereby depositing energy in the topology of the coro- 
nal magnetic field. Photospheric motions can also initiate coro- 
nal magnetic field disruption. Surface motions play an important 
role in forming Type B filament that are located between young 
and old dipoles and have a long, stable structure. This class of 
filaments (Type B) requires surface motions to gradually recon- 
figure preexisting coronal fields. 

In a previous study (Paper I), we removed all the large-scale 
flows in order to focus on smaller scale flow, such as mesogran- 
ulation and supergranulation. In this paper we have retained the 
large-scale flows in order to study their influence on the trigger- 
ing of a filament eruption. 

The three different methods used to estimate horizontal 
flows, while exhibiting small differences, provided a consistent 
picture of the general trends of the flow patterns. The LCT-B 
method gave smaller velocity amplitudes but the flow directions 
agreed quite well with the other methods. The smaller ampli- 
tude occurs because of the small number of the small isolated 
magnetic elements and the smoothing effect of the spatial win- 
dow used by this method. This leads to smoothing the results 
of the LCT-B method by approximately a factor of two. This 
agrees with the previous results of Schuck (2006), who shows 
that this method is not suitable for determining accurately esti- 
mating of the magnetic footpoint velocities. Correlation coeffi- 
cients comparing velocities from the manu-B and LCT-Doppler 
methods are positive and significant. General trends are very 
similar; however, we see many local discrepancies in the mea- 
sured flow fields due to the noise. This implies that magnetic 
elements (detected by the LCT-B and manu-B methods) do not 
necessary follow the plasma flow (detected by the LCT-Doppler 
method). 

The filament eruption started at about 16:30 UT at the lat- 
itude around -25° where the measurements of the horizontal 
flows based on Dopplergram tracking show a modification of the 
slope in the differential rotation of the plasma. This behaviour 
is not observed in the curves obtained by tracking the longitu- 
dinal magnetic field; both methods show a continuous slope in 
the differential rotation in the same place. This seems to be a 
consequence of the presence of a north-south stream along the 
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Fig. 8. A 3-hour average of the flow field in the immediate vicinity of the starting point (I — 56°, b = -26°) before (left) and after 
(right) the filament eruption. The filament observed by ISOON on 2004 October 7 at 15.20 UT is superimposed. The dashed boxes 
denote areas used for the zonal shear calculation. 




Fig. 10. Horizontal flow evolution over three days 6 {upper left), 7 (upper right), and 8 (bottom left) October, 2004. Bottom right - 
Mean zonal velocities, computed between -35° and +20° in longitude and -30° to 0° in latitude. Dashes (and triangles) represent 
mean zonal velocity for the LCT applied on Doppler on 2004 Oct 6. Solid lines (and squares) represent mean zonal velocity for the 
LCT applied on Doppler on 2004 Oct 7. Dash dots (and diamonds) represent mean zonal velocity for the LCT applied on Doppler 
on 2004 Oct 8. 
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filament position, which is easily measured by tracing Doppler 
structures and is only slightly visible in maps obtained by track- 
ing of magnetic elements. 

The observed north-south stream has an amplitude of 30- 
40 m s _1 . In the sequence of Ha image that record the filament's 
evolution, the part of the filament, which is in the north-south 
stream, is rotated in a direction compatible with the flow direc- 
tion of the stream. This behaviour suggests that the foot-points 
of the filament are carried by the surface flows. The influence 
of the stream is strengthened by differential rotation. We should 
keep in mind that the filament extends from -5° to -30° in lat- 
itude and that the northern part of the filament is subjected to a 
larger rotation than the southern part. 

The north-south stream, along with contribution from differ- 
ential rotation causes the stretching of the coronal magnetic field 
in the filament and therefore contributes to destabilizing the fila- 
ment. The topology of the north-south stream changed after the 
filament eruption, nearly vanishing. 

We have measured an increase in the zonal shear at the site 
where the filaments eruption begins, before and its sudden de- 
crease after the eruption. This result suggests that the shear in 
the zonal component of the flow field is the most important com- 
ponent of the surface flow affecting the stability of the coronal 
magnetic field, and it can lead to its eruption, which in turn can 
drive active phenomena such as ribbon flares and CMEs. This 
evolution of the shear in the flow field is probably related to the 
re-orientation by 70° (or 110°) of the transverse field after the 
eruption seen in the daily vector magnetograms obtained with 
THEMIS (Paper I)). 

All of the features observed in the topology of the horizontal 
velocity fields at the starting-point site could contribute to desta- 
bilizing the filament, resulting in its eruption. The present study 
has only examined the flows in the vincinity of a single filament. 
From our data, we propose that the stability and evolution of the 
filament are influenced by surface flows that carry the footpoints 
of the filament. In addition, Dudik et al. (2007) constructed a 
linear magnetohydrostatic model of the northern part of the fila- 
ment. His models show that the shape of the dipped field lines of 
the central part of the filament footpoints closely resemble the 
shape of the underlying, nearby polarities. This suggests a re- 
connection could be taking place between the flux of the incom- 
ing parasitic polarity and the "native" flux of the weak polarities 
dominant in that part of the filament channel. 

Filaments, or prominences, are important complex structures 
of the solar atmosphere. Several mechanisms are probably in- 
volved in the filament eruptions: the action of surface motions 
to create or increase the helicity of the flux rope (van Driel- 
Gesztelyi 2005, Romano 2005), reconnecting field lines in the 
corona (Mackay and van Ballegooijen, 2006a), the chilarity evo- 
lution of the barbs (Su et al. 2005), and oscillations of the fila- 
ment (Pouget 2006), etc. The coronal magnetic field is generally 
thought to be anchored in the photosphere, and flux transport on 
the solar surface (Wang et al. 1989) is the natural mechanism 
to explain the evolution of filament. Recent models of the large- 
scale coronal structure (Mackay and van Ballegooijen, 2006a) 
consider the action of the large-scale surface motions, such as 
differential rotation, meridional flow, and surface diffusion (su- 
pergranular). Recent analysis of the near-surface flows computed 
from Doppler imaging, provided by the MDI/SOHO instrument, 
reveals the shearing flow aligned with the neutral line (Hindman 
et al. 2006). Our present observation indicates that large-scale 
surface flows are structured (not uniform), showing areas of di- 
vergence or stream flows that should be taken into account in the 
numerical simulations. 



A better understanding of the mechanisms that lead to fil- 
ament eruptions requires simultaneous multi-wavelength and 
multi-spatial resolution observations (both high resolution of the 
filament and low resolution of the full sun) over a wide range of 
latitudes. Indeed, our previous works showed that different phe- 
nomena are observed at high resolution, such as magnetic re- 
connection close to the starting location of the filament eruption 
(Paper I). In this paper, observing a larger area at lower reso- 
lution, we showed that, at the same location where the filament 
first begins to erupt, there is a steep gradient in differential rota- 
tion, a north-south stream, and a shear in the zonal component. 

The next step in our study of the filament eruptions will be 
to examine the evolution of the extrapolated coronal field from 
photospheric longitudinal magnetograms to determine whether 
the effects of surface motions on coronal fields play a critical 
role in causing filament eruptions. 
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